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Abstract 

Using the theorem of residues Chiarella and Reichel derived a series 
that can be represented in terms of the complex error function (CEF). 
Here we show a simple derivation of this CEF series by Fourier expan- 
sion of the exponential function exp (— r 2 /4) . Such approach explains 
the existence of the lower bound for the input parameter y = Im[z] 
restricting the application of the CEF approximation. An algorithm 
resolving this problem for accelerated computation of the CEF with 
sustained high accuracy is proposed. 

Keywords: complex error function; Faddeeva function; Voigt func- 
tion; plasma dispersion function; complex probability function; com- 
plimentary error function; spectral line broadening 

1 Methodology and algorithm 

The complex error function (CEF), also known as the Faddeeva function, is 
given by 

w{z) = K{x,y) +iL (x,y) , 
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where z = x + iy and its real and imaginary parts are 



00 2 

K (x, y) = —= / exp I I exp (—yr) cos (xr) dr, y > 

V t J V 4 / 





and 



00 2 

£ 0, y) = ^= J exp ex P ( _ 2/ r ) sin ( a;r ) rfr ' 2/ > °> 


respectively. The real part of the CEF is known as the Voigt function. Com- 
bining the real and imaginary parts together yields 

00 2 

w (z) = —j= j exp ^— — ^ exp (—yr) exp (ixr) dr. (1) 


Let us show how the Chiarella and Reichel series approximation can be 
derived in form of the CEF [H [2] by Fourier expansion of the exponential 
function method. Consider the exponential function approximation that can 
be obtained either by Fourier expansion [31 H] or, equivalently, by Poisson 
summation formula [5] (see also work [6] for the Poisson summation formula 
application) 

/ t 2 \ a f n-K \ 

exp ( -— 1 « -— + a n cos I — r , -r m < r sC r m , (2) 

4 / ^ „= V Tm / 



where the Fourier expansion coefficients are 

2^ / n 2 7r 2 
a n « exp I — 

r m and iV are some parameters that may be chosen, for example as 12 and 
23, respectively. 

The approximation (j2J) is valid only within domain — r m ^ r ^ r m . Specif- 
ically, while the left side of equation (J2]) is non-periodic, its right side is pe- 
riodic (see Fig. 1 in Ref. [3]). Consequently, in order to obtain a higher 
accuracy we have to restrict the integration within the domain r £ [0, r m ] 
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after substitution of approximation §Z§ into integral (pQ) (see Ref. jl] for 
details), i.e.: 



w (z 



1 rn 

a 



2^ 



J exp (—yr) exp (ixr) dr 



o 

N 



(3) 



H — / a n cos( — r ] exp (— yr) exp (ixr) dr. 



i _ gi (ri7r+r m z) ^ — g' (-™ 7r + r m z ) ' 

/ , a„r m ; | — a - 

* — ' v nir + T m z mT — T m z 



The integral terms taken analytically in this approximation results to [I] 
r n 

that after trivial rearrangements using e ±tnn = (—1)™ can be simplified as 

Tm-Z V 71 " i ™ 7T - ^ 

It should be noted that the approximation (j3J) is valid only for positive 
value of y. However, due to symmetric properties of the real and imaginary 
parts of the CEF, its application can be easily extended for negative value 
of y as well (see for example Refs. [HE])- 

Since the CEF (Q]) contains the multiplier exp (—yt) that at relatively 
large y = lm[z] effectively damps the integrand to zero as r increases, we can 
assume that an extended integration will cause just a negligible error in com- 
putation. Based on this assumption we will try to extend the upper limit in 
integration to infinity. Thus the integrals above can be further approximated 
as 

T m OO 

J exp (—yr) exp (ixr) dr m J exp (—yr) exp (ixr) dr 



and 



; a n cos 

o 



T, 



oo 

' \ f f n7T \ 

-r exp (—yr) exp (ixr) dr ss a n cos I — r exp (—yr) exp (ixr) dr. 

: J J \T~rn J 
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Consequently, the series approximation can be rewritten in form 
2A. 



w (z) ~ j= I exp (—yr) exp (zxr) dr 



o 

oo 



1 N 7 (mx \ 

H — = > / a n cos I — r ] exp (—yr) exp (ixr) aV 



(5) 



Each integral term in approximation (jSJ) can be found analytically. This 
leads to 



( v . a . ,r^(a: + ij/) I 
w U ~ —i-—=r- — + z = > a r 



or, after slight rearrangement, to 
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2ir m zV^^ — . (6) 



r m z ' n-"K" - t*z* 

n=l 



Defining a small parameter h = vr/r m , the approximation can be 
expressed in a more traditional form 



W 2 ~ « 2 > -^T7> 9' 

ix z 7T z — ' n z h z — z z 

11=1 



This is the Chiarella and Reichel series approximation expressed in terms 
of the CEF according to the literature [HE]- It is interesting to note that 
this approximation was implicitly rediscovered afterwards. In particular, 
substituting it into following identity [TT] 



w (z) = e ^erfc (—iz) -v^ erfc (z) = e z2 w (iz) 

yields an approximation of the complimentary error function [T21 [T3] (see 
also equation 2.9 in the Ref. [2]): 

hze -z 2 ( I " e -nW 



7r \ z 2 ' n 2 h 2 + z 2 

n=l 



Since erfc (z) is directly proportional to w (iz) where iz = —y + ix, we can 
similarly explain why the this complimentary error function approximation 
works only at larger value of the input parameter x. 



4 



Table 1: The generated numbers corresponding to the real part of the com- 
plex error function (the Voigt function). 



X 


y 


Approximation Q 


Approximation ([5]) 


Algorithm 680 


0.01 


0.01 


9.887176929549550E- 


•1 


4.196286232960261E0 


9.887176929549547E-1 


0.1 


0.1 


8.884785624756435E- 


•1 


7.590865094856971E 1 


8.884785624756436E-1 


0.5 


0.5 


5.331567079121748E- 


•1 


5^3J516264696^^ 


5.331567079121750E-1 


1 


1 


3.047442052569125E- 


1 


Xu474420518lll^9E^ 


3.047442052569128E-1 


2.5 


2.5 


1.167371250446503E- 


■1 


1.167371250446503E-1 


1.167371250446503E-1 


5 


5 


5.696543988817698E- 


•2 


5.696543988817699E-2 


5.696543988817697E-2 


7.5 


7.5 


3.777752935845998E- 


•2 


3.777752935845999E-2 


3.777752935846000E-2 


10 


10 


2.827946745423245E- 


•2 


2.827946745423246E-2 


2.827946745423246E-2 


12.5 


12.5 


2.260351678541391E- 


•2 


2.260351678541392E-2 


2.260351678541391E-2 


15 


15 


1.882714532513676E- 


•2 


1.882714532513675E-2 


1.882714532513676E-2 



Table 2: The generated numbers corresponding to the imaginary part of the 
complex error function. 



X 


y 


Approximation ([J} 


Approximation © 


Algorithm 680 


0.01 


0.01 


1.108529605747765E- 


•2 


4.137187541585456E0 


1.108529605747726E-2 


0.1 


0.1 


9.433165105728508E- 


■2 


2.042540773410453E 1 


9.433165105728510E-2 


0.5 


0.5 


2.304882313844584E- 


1 


2^3J4774733^%7^E>1 


2.304882313844584E-1 


1 


1 


2.082189382028317E- 


■1 


1^082^^938202^^ 


2.082189382028316E-1 


2.5 


2.5 


1.079085859964814E- 


1 


1.079085859964814E-1 


1.079085859964814E-1 


5 


5 


5.583874277539103E- 


•2 


5.583874277539103E-2 


5.583874277539103E-2 


7.5 


7.5 


3.744329372959511E- 


•2 


3.744329372959512E-2 


3.744329372959514E-2 


10 


10 


2.813843327633689E- 


•2 


2.813843327633690E-2 


2.813843327633690E-2 


12.5 


12.5 


2.253130329137736E- 


■2 


2.253130329137737E-2 


2.253130329137736E-2 


15 


15 


1.878535427799564E- 


■2 


1.878535427799565E-2 


1.878535427799565E-2 



Tables [T] and |2] show the numbers generated by approximations (j4j), (J6]) 
and highly accurate well-known Algorithm 680 [91 [10] , given for comparison. 
The numbers that are failed in computation are struck out, the numbers with 
insufficient accuracies are underwaved and the numbers with relatively high 
accuracies are underlined. All other numbers are highly accurate. As the 
accuracy of the approximation ()6]) becomes high at y > 1, it is reasonable to 
apply it at larger values of the input parameter y in order to accelerate the 
computation. 

At smaller value of parameter y the decay of damping function exp (—yr) 
occurs at a lower rate. Consequently, due to periodicity of the function shown 
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on the right side of approximation (J2J), a smaller decay increases the error in 
calculation. This apparently explains why the accuracy of the approximation 
dOJ) drastically decreases as y decreases. The only way to resolve it is to 
increase r m (note that the period is 2r m ) and integer N determining the 
number of summation terms. Technically, this makes its practical application 
very inconvenient for the full coverage of the required range y > 10 -4 in the 
radiative transfer applications |14j . 

In order to resolve this problem, we can represent the approximation (J3J) 
as follows 



u> (z) 



- 2r„ 



-2t„ 



JV 



(-l) r 



(7) 



the common part the refining part 

There are two distinctive parts. The first part is nothing but the Chiarella 
and Reichel series approximation of the CEF (JHJ). Therefore it is the common 
part. Since the accuracy of approximation (J4]) is essentially higher than (jHJ), 
the second part can be called the refining part. Evidently, the Chiarella and 
Reichel series approximation of the CEF is embedded into approximation 



The algorithm is simple and implemented by turning on or off the refining 
part. Initially the program verifies whether or not y is smaller than some 
given value, say y < 1. If yes, the refining part is turned on and both 
parts are involved in computation. Otherwise, it is turned off and only the 
common part, i.e. the Chiarella and Reichel series approximation of the CEF 
alone, remains active in the program. Such technique enables accelerated 
computation as the refining part may not be required when the value of 
y is sufficiently large. The described algorithmic implementation of CEF 
approximation in the form of (JTj) provides more rapid computation with 
sustained high accuracy. 



2 Conclusion 

The Chiarella and Reichel series approximation of the CEF is derived by 
Fourier expansion of the exponential function exp (— r 2 /4). This derivation 
methodology helps understand why application of this CEF approximation 
is restricted at smaller value of the input parameter y = lm[z\. An algorithm 
resolving this problem for accelerated computation with sustained high ac- 
curacy is proposed. 
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